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We describe the equation of motion of two charged spherical shells with tangential 
pressure in the field of a central Reissner-Nordstrom (RN) source. We solve the problem 
of determining the motion of the two shells after the intersection by solving the related 
Einstein-Maxwell equations and by requiring a physical continuity condition on the shells 
velocities. 

We consider also four applications: post-Newtonian and ultra-relativistic approxima- 
tions, a test-shell case, and the ejection mechanism of one shell. 

This work is a direct generalization of Barkov-Bclinski-Bisnovati-Kogan paper. 

Keywords: Classical gravity; Exact solutions. 



1. Introduction 

The mathematical model that we analyze in this paper describes the dynamic evo- 
lution of two spherical shells of charged matter which freely move outside the field of 
a central Reissner-Nordstrom (RN) source. Microscopically these shells are assumed 
to be composed by charged particles which move on elliptical orbits with a collective 
variable radius. The angular motion, distributed uniformly and isotropically on the 
shell surfaces, is mathematically described by a tangential-pressure term in the en- 
ergy momentum tensor of the Einstein equations. The definition of the shell implies 
that all the particles have the same following three ratios: energy/mass, angular 
momentum/mass, and charge/mass. Indeed, since the equations of motion for any 
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singled-out particle "a" are 
dt a 1 



ds -m a c 2 g tt (r a ) 



(E a + e a A (r a )) (1) 



dr a V 1 ,„ 1 \ / ll 1 \ 1 



-{E a + e a A (r a )) 



ds J m 2 a c A \-gu(r a )grr(r a ) J \m 2 a c 2 r 2 J g rr (r a ) 

(2) 



ddaV ll 1 



ds J m 2 c 2 r 4 m 2 c 2 r 4 sin 2 ( 
dip a k a 1 



(3) 
(4) 



ds m a c r 2 sin 2 a 

(gtt and g rr are the components of a spherical symmetric metric and A is the 
electric potential; k a and Z Q are arbitrary constants), it is easy to see that the radial 
motion for all particles is the same if 

— — = const, —2- = const, = const, Va, (5) 

m a m a m a 

where each const, does not depend on the index a. Therefore, if at the beginning the 
particles are on the same radius r a = Rq, then the shell will evolve "coherently", 
i.e. all particles will evolve with the same radius. 

Now the problem we are interested in is to find the exchange of energy between 
the two shells after the intersection. Indeed the motion of the shells before and 
after the crossing can be easily deduced from the equation of motion for just one 
shell, which equation has been found many years ago by Chased with a geometrical 
method first used by IsraeP. Instead, the intersection problem was considered first 
by Langer-Eiol they applied it to the particular case of electrically neutral and 
pressureless shells (neutral dust), then Barkov et al W considered the more general 
case of shells intersections when the shells have also tangential pressure (so Langer- 
Eid's results follows from Barkov et al. results as particular case). Now we generalize 
this problem for the shells with tangential pressure and also electric charge. 

What we achieve in the present paper is the determination of the constant pa- 
rameters after the intersection knowing just the parameters before the intersection. 
Actually the unknown parameter is only one, mix, which is the Schwarzschild mass 
parameter measured by an observer between the shells after the intersection. This 
parameter is strictly related to the energy transfer which takes place in the crossing, 
and it is found imposing a proper continuity condition on the shells velocities. 

In the model we assume that there are no other interactions between the two 
shells apart the gravitational and electrostatic ones. In particular the shells, during 
the intersection, are assumed to be "transparent" each other (i.e. no scattering 
processes) . 

The paper is divided as follows: in Sec. 2 we preliminarily discuss the one-shell 
case; in Sec. 3, which is the central part of this article, we find the unknown param- 
eter 77121; then, Sees. 4-7 are devoted to some applications: post-Newtonian approx- 
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imation, zero effective masses case (i.e. ultra-relativistic case), test-shell case, and 
finally the ejection mechanism. 

In this paper we deal only with the mathematical aspects of the problem; some 
astrophysical applications of charged shells in the field of a RN black hole have been 
considered in Ref. [3]. 



2. A gravitating charged shell with tangential pressure 

The motion of a thin charged dust-shell with a central RN singularity was firstly 
studied by De La Cruz and Israel^, while the case with tangential pressure was 
achieved by Chase^ in 1970. All these authors used the extrinsic curvature tensor 
and the Gauss-Codazzi equations. However we followed a different way, indeed the 
same solution can be found also by using 5 and distributions and then by di- 
rect integration of the Einstein-Maxwell equations (see Ref. [6] and the appendix in 
Ref. [7]). This method has the advantage of a clearer physical interpretation, and 
it is also straightforward in the calculations; however in the following we will give 
only the main passages. 

Let there be a central body of mass mm and charge ej„ and let a spherical 
massive charged shell with charge e move outside this body. It is clear in advance 
that the field internal to the shell will be RN, while externally we will have again a 
RN metric but with different mass and charge parameters m ou t and e ou t = e,-„ + e. 
Using the coordinates x° = ct and r, which are continuous when passing through 
the shell, we can write the intervals inside, outside, and on the shell as 



- (ds)l = -e T ^f m (r)c 2 dt 2 + fr^r)dr 2 + rW (6) 

- (ds) 2 out = -fout(r)c 2 dt 2 + f^ t {r)dr 2 + rW (7) 
-(d S ) 2 on = -c 2 dT 2 +r (T) 2 dn 2 (8) 



where we denoted 



dW = dQ z + sin' 



and 



in out . G( 

Jin = J- - ^ 5 ' 4 9~ i Jotlt = l — * n 1 4 — 5 ■ l y j 



In the interval (JHJ>, t is the proper time of the shell. The "dilaton" factor e T ^ in 
([6]) is required to ensure the continuity of the time coordinate t through the shell. 
If the equation of motion for the shell is 

r = Ro(t), (10) 
then joining the angular part of the three intervals (JHJ-fH]), one has 

r (r) = Ro[t(r)}, (H) 
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where the function t(r) describes the relationship between the global time and the 
proper time of the shell. Joining the radial-time parts of the intervals ©-([7]) on the 
shell requires that the following relations hold: 

/,Jr ,(§)V<'>-/-Hr„)(^) J = l, (12) 



/~M(£)'-i=iM (£)'_,. M 

If the equation of motion for the shell — i.e. the function ro(r) — is known, then 
the function t(r) follows from (Til?)) and consequently T{t) can be deduced by (fi"2"j) . 
Thus the problem consist only in determining ro(r), which can be done by direct 
integration of the Einstein-Maxwell equations 

R k i - \R9 k i = ^t? 

(14) 

with the energy-momentum tensor given by: 

T? = 6 Ul u k + (S 2 6 k + 5!S k 3 )p + k (15) 
T (ei) k = ±-(F a F kl - U k F lm F lm ) . (16) 

Here on we employ the following notations: 
— ds 2 — gikdx l dx k , gik has signature (— , +, +, +) 

x k = (ct,r,0,tp) i,j,k... = 0,1,2,3 

p = p(Ro) = pe = Pip ^tangential pressure (p r = 0) 

Fik = — 

The above equations are to be solved for the metric 

- ds 2 = g o(t, r)c 2 dt 2 + g n (t, r)dr 2 + r 2 dn 2 , (17) 

and for the potential 

A = A (t,r), A 1 =A 2 =A 3 = Q, (18) 

As follows from the Landau-Lifshitz approach [S] (see Ref.[B]) the energy distri- 
bution of the shell is 

M{t)c 2 5[r - Ro(t)] 



47rr 2 M°V-5oo5ii 

while its charge density is 

ceS[r-Ro(t)] 
P 47rr 2 M°^/-5oo3ii ' 



(19) 



(20) 



January 26, 2009 13:19 WSPC/INSTRUCTION 

Charged'Shell'Intersections 



FILE 



Intersections of self- gravitating charged shells 5 

where 5 is the standard 5-function. In the absence of tangential pressure p, the 
quantity M in Eqn. (fT9| would be a constant, but in presence of pressure, Mc 2 
includes the rest energy along with the energy (in the radially comoving frame) of 
the tangential motions of the particles that produce this pressure. 

It can be checked that the Einstein part of (|T4]) actually lead to the solution 
©-([8]) with, in addition, the "joint condition" 




v ' /m(ro)+ (^> +^-'( r °)+(^ = 2 imou *,- mm) ~ e2 i?z m > (2i) 

where we denoted 



while 



cdr J /j,(t) ^(r)c 2 ro 

M (t) = M[t(r)}, (22) 
niout - rriin — E I c 2 (23) 



is a constant which can be interpreted as the total amount of energy of the shell. 
Then, from the Maxwell side of (|14p the only non- vanishing component of the electric 
field is 



Fqi = _V_J[ooffii {e . n + ed[r _ Rom (24) 

(9(x) is the standard step function). Finally, the equations T£ k = can be reduced 
to the only one relation: 

p _ dM c 2 5[r-R (t)] 
dt 87rruV-foo5ii 

We will not treat here the steady case (i.e. ro = const) which should be treated 
separately; thus in the following we will assume always ro ^ const. . 

The joint condition l|2ip can be written in several different forms: two of them, 
which will be useful in the following, are 

/, , s , f dr a \ 2 {m out -m ln ) G^ 2 (r) - e 2 - 2ee in 
V /fc(n,)+ W = W) + (26) 



and 



/ A . f dr \ (m ov 



Gti 2 (T) + e 2 + 2ee v 



(27) 



(jl(t) 2^(r)c 2 r 

As in Ref. [BJ all the radicals encountered here are taken positive, since for astro- 
physical considerations only these cases are meaningful. To proceed further, we must 
specify the equation of state, i.e. the function (a(t). Here we consider a particle-made 
shell, therefore the sum of kinetic and rest energy of all the particles is 




Mc 2 = > K^l + -Ts - (28) 
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where p a is the tangential momentum of each particle (the electric interaction be- 
tween the particles is already taken into account by the self-energy term of, e.g., 
(|26"]) . thus one has not to include it in M too). From the definition of the shell (see 
Introduction) it follows: 

_P|_ ^ l l = const 

ml mlRl Rl ' ^ 
the square root in (|28p does not depend on the index a; then defining 

^m a c 2 =mc 2 , ^ \lg\ = L, 

a a 

formula (|28p can be re-written (remembering definition (f^|) too) as 




^) = \l m2 + ^ry (3°) 

Thus, now, one can determine the function tq(t) from equation (|21|) [or from one of 
the equivalent forms (|26 | - ([27|) ] if the initial radius of the shell and the six arbitrary 
constants m m , m out , m, e^, e and L are specified. Accordingly with (|19|. (|25| . 
(|22|) and (f30|). the equation of state that relates the shell energy density e to the 
tangential pressure p is 

L " - (l + — - 

?o V m 2 c : x l(1 



as in the uncharged case, i.e. the presence of the charges do not modify the relation 
between energy density and pressure (indeed the presence of the charge is hidden 
in the equation of motion). Note that when the shell expands to infinity (i?o — > oo) 
the angular momentum becomes irrelevant and the equation of state tends to the 
dust case p << e. 



3. The shells intersection 

Let us now consider the case of two shells which move in the field of a central charged 
mass. The generalization from the previous (single-shell) case is straightforward if 
the shells do not intersect: indeed the outer shell do not affect the motion of the inner 
one, while the inner one appears from outside just as a RN metric. Therefore the 
principal aim of this section is to consider the intersection eventuality and to predict 
the motion of the two shells after the crossing, having specified the initial conditions 
before the crossing. After the intersection one has a new unknown constant that 
has to be found by imposing opportune joining conditions as now we are going to 
explain (the analysis follows step by step the Ref{6js one). 

Let us previously analyze the space-time in the (t, r) coordinates (which are 
continuous through the shells). We define the point O = (t*,r*) as the intersection 
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Fig. 1. The four region in which it is divided the spacctime; the two lines represent the trajectories 
of shcll-1 and shell-2. 

point; then the space-time is divided in four regions (see Fig.l): 

COB (r > i?i,r > R 2 ), 
CO A (R 1<r < R 2 ), 

AOD (r < Ri,r < R 2 ), [ ' 

BOD (R 2 <r< R\). 

Correspondingly to these regions we have the metric in form (| 13|) but with different 
coefficients <?oo and gn- 



(COB) 

9oo 


= -fout(r) , 


(COB) 

an 


= fou\(r) 


(33) 


JCOA) 
#00 


= -e Tl W/ 12 (r), 


JCOA) 

an 


= /r» 


(34) 


„{AOD) 

aoo 


= -e To «/ i „(r), 


AAOD) 

an 


= /*(»•) 


(35) 


(BOD) 

aoo 


= -e T ^f 21 (r), 


JBOD) 

an 


= / 2 i 1 W 


(36) 



The dilaton factor T{t) allows to cover all the space-time with only one ^-coordinate; 
here, fi n and f out are the same as those in ^ while f\ 2 and f 2 \ are given by similar 
expressions: 

2Gmi2 G(e m + d) 2 
•H2 - 1 ^ 1 ^2 ( 6 '> 

2Gm 21 G{e m + e 2 ) 2 

/21 = 1 2 1 \ 6ti > 

c r c r 

As we said, the parameters n%i n , m\ 2 , m out , em, ei, e 2 are assumed to be specified 
at the beginning, while m 2 \ is the actual unknown constant which has yet to be 
determined from the joining conditions on (t», r*). 
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Before the intersection 



Let us write the equation of motion for the two shells before the intersection (shell- 1 
inner and shell- 2 outer). This can be made easily adapting the (|2"7|) and (f2"o) to the 
present case: 



, dn \ (mi 2 -m in ) GM\ + ef + 2e in ei 
/u(ri) + -7- = 77 TTT^I ( 39 ) 



cdri/ Mi 2M 1 c 2 r 1 

for shell 1, while for shell 2 



dr 2 \ 2 (m 12 -m ln ) GM% — e\ — 2{e in + ei)e 2 



' ,4o) 



with 



M i = y m i + ^2' M2 = y m * + ^f- (4i) 

Here, ri and r 2 are the proper times of the first and second shells respectively, 
while ri(ri) = Ri[t(ri)] and r 2 (r 2 ) = i^tfa)]. Now we have to impose the joining 
conditions for the intervals on both the shells. For the first shell (on curve AO) one 
has: 

^-w^)'-^"'^)'- 1 ' (43) 

while for the second shell: 

'-<■»> (k)'-^ (as)'" 1 (44) 



If all free parameters and initial data to Eqs. (|39|) - (j41|) were specified and if the 
functions fi(ri) and r 2 (r 2 ) were derived, then their substitution in f|42[) - (|45 () gives 
the functions Ti(t), r 2 (i) and Ti(t), To(t), which is enough for determining the 
motion of the shells before the intersection. Therefore the intersection point (t*,r*) 
can be found by solving the system 

r* = r 2 (r 2 (Q) , {W) 



which we assume that has a solution. 
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After the intersection 

The equation of motion for the shells after the intersection time i* can be con- 
structed in the same way again by turning to Eqns. (126|) and (|27p . and introducing 
the new parameter m 2 \ which characterize the "Schwarschild mass" seen by an ob- 
server in the region BOD. We use Eq. (f2"6"|) for (now outer) shell 1 and Eq. ([2T|) for 
(now inner) shell 2: 



, . drx\ (m out -m 2 i) , GMl - e\ - 2ei(e in + e 2 ) 

/2i(ri) + (^r) = — m[ — + 2ma ' (47) 



/, , v , / dr 2 \ 2 {m 2 i-m in ) GMf + e\ + 2e 2 e in 

V /2l(r2)+ feJ = M 2 2M^ ■ (48) 

Naturally, M\{r-\) and M 2 (r 2 ) are given by the same expression of (|4ip but now 
they have to be calculated on ri(ri) and r 2 (r 2 ) after the intersection. 
Joining the intervals on the first shell (on curve OB) yields 

'^(k)'-*©'- 1 (49) 

e Tm Mri) (^)\ J - Hri) (iiL) 2 = 1 . (50) 



k dri J \ cdri 

Then, joining the second shell (on curve OB) we obtain: 

Since the initial data to Eqs. (|4"Tj) and (|48[) have already been specified (from the 
previous evolution), then the evolution of the shells after the intersection would 
be determined from Eqs. ()47[) -(|52 p if parameter m 2 \ were known. Thus we need an 
additional physical condition from which we could determine m 2 \. 

This condition follows from the fact that the Christoffel symbols (i.e. the accel- 
erations) of the shells have only finite discontinuities (finite jumps), therefore the 
relative velocity of the shells must remain continuous through the crossing point. 

In the presence of two shells, we can construct one more invariant than in the 
single shell case (where only UiU 1 = — 1 was possible): the scalar product between the 
two 4- velocities of the shells. We can also avoid to apply the parallel transport if we 
evaluate the 4- velocities on the intersection point (i*, r*). The continuity condition 
can be found imposing that the scalar product has to have the same value when 
evaluated in both the two limits t — > tl and t — > tf. 
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Determination of Q. 

Let us start determining the quantity 

n _ { (COA) o „0 , JCOA) 1 i ; ^ 
^ = ISOO u AO u CO^~9ll u AO u COSt=t* ,r=r x =r 2 =r*, (oo) 

which is the scalar product of the two 4- velocities evaluated in the intersection point 
from the region AOC (along the curves AO and CO). Written explicitly, the unit 
tangent vector to trajectory AO is 

u AO = ( u AO' U AO' U AO> u Ao) 

dt dri \ . 
-j-.-Ao.O , (54) 

Ai cdr i Jt<u 

while for the trajectory CO we have 



u ho — ( u co> u co> u co> u co) 

dt dro \ , 

—1,0,0 . (55) 
dr 2 cdr 2 J t <t, 

The fact that these are actually unit vectors follows from the joining equations 
62D and HU). 

The components of the vector (|54|) can be easily expressed from Eqs.(|39|) and 
(|42| as 

e" Tl(t)/2 f GM 1 2 (r 1 ) + e j+2e 1 e m ^ 

mi2 - m» n ^-o (56) 



dnA< t Mt^hiin) V 2c 2 n 



CC?Ti 



t<t» 



GM 1 2 (r 1 ) + e 2 + 2e 1 e m . 
mi2 - m in —J - Mf (n)/i 2 (ri) 



M 1 (r 1 )/ 12 (r 1 ) V V 2c 2 n 

(57) 

where 

5 1 = s S n(^) . (58) 

Analogously, for the components of vector (|55p , we obtain the following expressions 
from Eqs.gQl) and 63) : 



dt\ e -ri(t)/a ^ GM|(r 2 ) - e 2 - 2e 2 (e m + ei) 



dr 2 J t < u M 2 (r 2 )/ 12 (r 2 ) V 2c 2 r 2 



(59) 
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dr 2 \ 
cdT 2 J t <u 



M 2 (r 2 )fi2{r 2 ) 



m out - mi2 + 



GM 2 2 {r 2 )-e 2 2 -2e 2 {e m + e 1 ) 
2c 2 r 2 



M|(r 2 )/ 12 (r 2 ) 
(60) 



«— - (asL- (6i) 

Thus, from the preceding results, we obtain: 

Q = M1M2/12 ' 

fY m m Gj\/ 1 2 +e;+2e iei „ \ / GM| -e| -2e 2 (e„, + ei ) \ 

re // GM 2 +ef+2eie in \ 2 , f2 , 

-oifay ( - m in 1 2c i rr J -Mf/12 

// , GM 2 2 -e|-2e 2 (e 1 „+ei)\ 2 , , r2 , \ 

y ("W - mi2 H 2 — 2 2 ?7 ^ -J -M|/i 2 j; 

(62) 

here and in the following we omit the coordinate dependence of f a , M a etc., im- 
plicitly assuming that they have to be evaluated on (t*,r«) where not differently 
indicated. 



Determination of Q' . 

It is possible to apply the same procedure to the region BOD (i.e. after the inter- 
section time), finding the quantity 

Ql = r (BOB) , (BOD) 1 1 , 

V — lf00 u OB u OD a OB a ODit=t* ,r=r 1= r 2 =r, ■ (Vd) 

Now the unit tangent vectors to trajectories OB and OD ar^f]: 



i bs — ( u ob> u ob> u ob> u ob) 
dt dr\ 



dr\ ' cg?ti 



,0,0 



(64) 



t>t* 



and 



a Obviously, when we say t > t*, we tacitly assume before a (possible) second intersection. 
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u OD — ( u ODi u ODi u ODi u Od) 

P-^Ao) ; (65) 
dr 2 cdr 2 ) t > u 

from the joining conditions (|50| and (|5Tjl it is possible to see that these are 
actually unit vectors. The components of u l OB can be deduced from Eqs.(j47| and 
(f50|) . while the components of u % OD from Eqs. pS)) and (|5T|) . Then, using the metric 
in the region BOD, it is possible to calculate the scalar product 



M1M2/21 

GMn+ei+2e 2 ei. 



2c 2 r„ 



m ou t - m 2 i H J — 2c 2 r, J ( m 2 i - m m - 

XI X' If . GM 1 2 -e 2 -2ei(e IIl +e 2 )\ 2 . f2 , 

~"i"2V y^out — mix H 1 — 1 - 2 -^ TF f -J -Mf/21 

™ GM 2 + e 2 +2e 2ei „ \ 2 M2f \ 

m 2 i — m in 2^7, ) ~ M 2 /21 r , 

(66) 

where S[ and 5' 2 have been defined as in (|58j) and (f61|) . but for t > t*. We introduced 
these symbols only for generality, but actually we are interested only in the case 
withEl 

6[=S U 5' 2 = 5 2 . (67) 
The necessary continuity requirement is thus 

Q = Q', (68) 
then, since r* is assumed to be known, this equation allows to find 77121. 

Physical meaning of Q and Q' . 

Using standard definition for the shell velocities before the intersection one has 

1 2 _ 9 , n° A) (n) (drx 



"i> 2 \ 2 _ g[i° A \r 2 ) f dr 2 



-g£ OA \r>)\cdt 



(70) 



b This is the only possible case if one excludes vi(t*) = V2(t*) = 0, because there are non discon- 
tinuities in the velocities. 
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and similarly for the velocities after the intersection, 

' ,\2 {BOD), x / , x 2 



g{i° D \r 2 ) ( dr 2 



Then it is easy to obtain from the definitions and (|63p. thajf) 



(71) 
(72) 



y/l ~ Vl/C^l - V$/C 



2 



t— 1„ ,ri — T2— 



and 



Q' = - . (74) 

Determination of P and f . 

First of all it is convenient to introduce new symbols to simplify the expressions of 
Q and Q'. With 

GM\ + e\ + 2eie in 
Qi = 



'12 



2c 2 r* 

_ GAf 2 2 -e|-2e 2 (e m + ei) 



2c 2 r 



and 



, GAf 1 2 -e 2 -2e 1 (e m + e 2 ) 

91 " 2c 2 r* 

, _ GM 2 + e 2 + 2e 2 e m 
12 — 



2c 2 r 

then Q and Q' can be re-written as 

Q = A/i A/2/12 ' 



• j(mi 2 - m in + gi) {m out - m, 12 + g 2 ) + 
-S1S2J (mi2 - m ln + qif - M 2 / i2 
\J (m out - m 12 + q 2 ) 2 - M|/i2 j 



(75) 



c It is also worth noting that ^/ Q 2 — 1/Q = —\v\/c — V2/c\/(l — v\V2/c 2 ), which is the relative 
velocity definition of two "particles" in relativistic mechanics. 
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and 



Q' 



M1M2/21 



(76) 



■Urriout ~ m 2 i + q[) (m 2 i - m in + q' 2 ) + 

-$'i&2\] («W - rn 2 i + q'x) 2 - Mffo 

\J (m 21 - m in + q 2 f - M$f 12 } , 

Now, in principle is possible to find 77721 by squaring and solving Q = Q' (which 
is a quartic equation). However the procedure is cumbersome and moreover it is 
not possible with Eq. ([68|) alone to determine the sign of the roots. Fortunately, 
as in the non-charged case, it is possible to follow another easier way. Indeed, it 
is possible to introduce two other invariants, say P and P', similar to Q and Q' , 
which are constructed using the scalar products of the 4- velocities of the two shell, 
but now taking the limit to (i»,r») from the AOD and COB regions respectively. 
More explicitly, we define 

p = t {AOD) , (AOD) i 1 , (m 

r — ti/00 a AO a OD * 9ll a AO a ODH=t* ,r=r 1= r 2 =r* , {II) 

and 

p/= r (cob) a (cob) x 1 , (m 

r — IffoO u CO U OB+9ll u CO U OBJt=f ,r=r 1 =r 2 =r t ■ (<o) 



Then, the same continuity requirement of Eq. (|68[) implies that it must hold also 
that 

Q = P , P = P' . (79) 

Following the same method used to find Q and Q', after some calculations, one 
arrives to 

P = 



M 1 M 2 fir, 



•j (77112 ~ 777i„ +Pl) (77721 - 777m + P2) + 

-8i5' 2 \J (77712 - m in +Pif - Ml f in 

\J (77721 - 777i„ + p 2 f - M|/j„ | 



(80) 
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and 



P' 



-l 

Af 1 M 2 / < „ 



■|(m ou t - m 2 i +p[) (rriout - m\i + p' 2 ) + 
-S'ifoyJ {m ou t - m 2 i +p'i) 2 - M\f out 
\j (rn out - mi2 + p' 2 ) 2 - M$f out j , 

where we have denoted 

GM\ -e\- 2 ei e in 



(81) 



Pi 

P2 



2c 2 r« 
GM 2 —ek — 2e 2 e in 



2c 2 ? 



and 



, _ GA/ 1 2 + e 2 + 2ei(e m + e 2 ) 
Pl ~ 2c 2 r» 

, _ GM 2 +e| + 2e 2 (e m + ei) 
P2 ~ 2c 2 r» 

Determination of m?i; the energy transfer 

Thus the complete set of continuity conditions at the point of intersection can be 
written as 

Q = Q\ Q = P, Q = P'. (82) 

It turns out that this three quartic equations for the unknown parameter m 2 i have 
only one common roo10 . It is possible to find the solution using hyperbolic functions. 
The final result is remarkably simple: 

ei e 2 GM1M2 

m 2 i = m in + m out - m 12 = = Q , (83) 

err* c z r» 

or equivalently, in terms of f 2 \. 

Hi = fin + font - fa + 2 — 2 Q ■ (84) 

c r* 

It can be easily seen from Eqn. ([83|) that the charge e; n of the central singularity 
does not affect the result (but it affects the equation of the motion of the shells and 



d In the pure gravitational case it is also possible to use just Eqn. Q = P' , and then choose the 
correct root solution by requiring that m\2 has to be positive, see Ref. [5] 
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thus Q). Formula (|83p solves the problem of determining the mass parameter m 2 i 
from the quantities specified at the evolutionary stage before intersection. It is then 
possible to determine the energy transfer between the shells. Indeed the energy of 
shell 1 and 2 before the intersection are, respectively 

Ei = (TO12 - m in )c 2 , E 2 = {m out - mi 2 )c 2 , (85) 

while, after the intersection 

E[ = (mout ~ m 2 i)c 2 , E! 2 = (m 2 i - m in )c 2 . (86) 

The conservation of total energy is automatically ensured by the above formulas, 
indeed 

Ei + E 2 = E[ + E' 2 . (87) 
Then it is natural to define the exchange energy as 

AE = E' 2 -E 2 =-(E[ -Ei) . (88) 
Then, from Eqn. (f83|) and the above definitions, it follows that 

AE = - ^MlQ . (89) 

It is also useful (especially for the Newtonian approximation) to use Eqn. (|T3")) and 
re-express AE as: 

Ap _ eie 2 GMiM 2 \ viv 2 /c 2 -l \ 

r* \ y/l-vi/fy/l-vi/c? j r=rt 

4. Post-Newtonian approximation 

For slow velocities of the shells it is interesting to consider the Post-Newtonian limit 
of Eqn.(|9T)l): 

^ E _ Gmim 2 - eigg | 

1 {Gmvm 2 , . . . , l9 Gm 2 L 2 GmiL 2 ,~\ ( \\ 

+ 7r 2 { — — [Mr* )-v 2 (n 2 + — V + — -r \ + - • ( 91 ) 

It is worth noting that only the zeroth order in 1 /c 2 changes with respect to the 
uncharged case (because of the Coulomb term — eie 2 /r*), while all the other orders 
remain unchanged, being of kinetic origin; mj and m 2 are the rest masses of the 
shells, indeed we have used for the masses Mi and M 2 the definitions (|4Tj) . 

It can be also useful to re-express all the quantities in a Newtonian language 
and consider only the zeroth order in 1/c 2 , e.g. we can expand the energy as 

E = mc 2 +£ + o f-M , (92) 
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where m and £ do not depend on c. Therefore, similarly, we can define at the first 
order in 1/c 2 

m 12 - m in = mi + -i, m out - m\ 2 = m 2 + -f-, (93) 
c c 

£ ! £' 
- m 2 i = mi + m 2i - m m = rai + — . (94) 

c c z 

Then it follows also that the energy conservation law takes the form 

£ i + £ 2 = Si + 4 , (95) 

and Eqn. ([55|) becomes 

£[ = £%- A£, £ 2 = £ 2 + A£, (96) 

where A£ = (Ai?) c ^ 00 . Thus from the above formulas and definitions it is clear 
that 

A£ = Gmim2 - 6162 . (97) 



5. Pressureless shells with zero effective masses (L± — L 2 — and 
Mi = M 2 = 0) 

It is interesting also to consider the case in which the motion of the particles of the 
shells is only radial (i.e. L\ = L 2 = 0) and the rest masses are negligible with respect 
to the kinetic energies and to the charges — indeed this is the case for two shells 
composed by (ultra)rela-tivistic electrons and positrons. In this case the effective 
masses can be replaced by 

Mi = M 2 = A , (98) 



where A is a parameter arbitrary small. From Eqn.(|89[). with Q expressed by formula 
([62"]) . it is easy to find that the energy transfer in this case is 

AE = - P -^ + -^H/ m - /i 2 )(/i2 - font) + o(X 2 ) , (99) 
r* 2G/i2 

having assumed that the shells have opposite-directed velocities, i.e. 

§ x 8 2 = -1 . (100) 
Otherwise, if the shells goes in the same direction, i.e. 

<Ji*a = 1 , (101) 

then Eqn. (|55|) becomes simply 

AE = -— + o{\ 2 ) ; (102) 

obviously the previous formulas make sense only if r* exists. We want to underline 
the presence of the term o(A 2 ), because, strictly speaking, a charge cannot have zero 
rest mass, therefore we are in the case of just small effective masses. As expected, in 
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the case of vanishing charges (ei = e 2 = 0), Eqn. (|102p gives zero at A = because 
this is the case of two photon-shells which go in the same direction and therefore 
cannot never intersect. 

6. The intersection of a test shell with a gravitating one 
One-shell case 

Let us consider firstly the case of a test shell on the RN field. This limit has the 
only aim to show that the shell's equation of motion (|26|) actually reduce to the 
simple test-particle case; the limit can be obtained by putting 

m — > Am , e — > Xm , L — > XL , (m out — m in )c 2 — » XE (103) 

with A — > 0. Then, considering also (f3"0")) . we find that Eqn. (f2"6")) becomes 



now, putting A = the self-energy term is killed; then re- writing Eqn. (1104p using 
the more familiar Schwarzschild time t and Eqn. ljlip . 



it is easy to recognize that Eqn. (|105p coincides with the first integral of motion 
of a test-charge particle on the Reinssner-Nordstrom background, where E is the 
conserved energy of the particle, m the rest mass, e the charge and L the angular 
momentum. 

Two-shell case, with one test-shell 

Now we can deal with the more interesting two-shell case, in which shell-2 is con- 
sidered "test" . To gain this limit we have to put 

m 2 — > Xm 2 , e 2 — > Xm 2 , L 2 — > XL 2 , (106) 

(m out - mi 2 )c 2 — ► Ai? 2 , (m 2 i - m in )c 2 — ► A£ 2 . 

Then, using Ean.(|83p with Q given by formula (j6"!?|) , one obtains 

AE = + -k- ■ 



2c 2 r, 



(107) 



(12 



+ A^) -M|/ 12 }. 
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Thus, only the self-energy terms of shell-2 are killed by A = 0. 

Now, it is worth noting the following fact: shell-1 does not have any discontinuity 
when it intersect the shell-2 (this is natural because shell-2 is "test" and does not 
affect the metric) , on the other hand shell-2 undergoes a discontinuity in the metric 
when it cross shell-1 and consequently it has an actual discontinuity in the velocity. 
It is easy to calculate this gap; indeed using the definition (|70|) of velocity v 2 [with 
the time given by the joint condition (|l5|) ]. with metric coefficient (fM)) . and 
with the help the first integral of motion (|40|) . one finds 

vl(r 2 ) = 1 - f out (r 2 ) f-^L- - e f/. + " m)2 ) 2 + o(A) , t < U , (108) 
\M 2 (r 2 ) M 2 (r 2 )r 2 J 

where we have used f\ 2 = f ou t + o(A); in the same way, using (|72l) . (|5T|) . (|36j). and 
(|48| . the velocity v 2 (after the intersection) is 

where E' 2 can be expressed in function of E 2 with the help of (|107j) . From the 
previous formulas it is clear that in general 

v' 2 (r*) - v 2 (r*) t^O . (110) 



7. Shell ejection 

The exchange in energy of the shells during the intersection makes possible that 
one initially bounded shell can acquire enough energy to escape to infinity. 

The shell ejection mechanism can take place also in the Newtonian regime. In 
this case, from Eqs. (f96"|) - (f97)) it results that 



, Gm 1 m 2 -e 1 e 2 , Gm 1 m 2 -e 1 e 2 

t 1 — tl ; , t 2 — t 2 H , (111J 

and then, after the first intersection 



ell CI Gm 1 m,2—e 1 e2 c 

L 2 — fc 2 —i — 



2 + (Gm 1 m 2 - e x e 2 ) (jr - ^rrj , 



(112) 



where we have denoted the radius of the first and second intersection with r' r and 
r'l respectively. In the following we will consider only the case 

Gm\m 2 — e\e 2 > , (H3) 

this is e.g. the case in which the two shells have opposite charges. Thus, also in the 
case £i,£ 2 < 0, if 

r'l > K » (H4) 
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and if the initial condition were in such a way that rj, is enough small and r" not 
too much close to r*, then it is possible to have £ 2 > 0, i.e. the ejection of the 
second shell. 

Let us now assume that r" > , and consider a "semi-relativistic" case in which 
at the first intersection we use the full relativistic formulas^, 

M 1 (r',)M 2 (r'„] 



E[ = Ei 



E — E 2 



H-Q) 



(115) 



MiK)JIJ 2 «) 



(-Q) - , 



while at the second intersection we use the Newtonian approximation, 



E'l 



rpl I Gra\rn 2 — e\e 2 
1 r'J 

£j _ \ M 1 (r' t )M 2 (r' t )(-Q)-e 1 e 2 _ Gm 1 m 2 -e 1 e 2 



J 2 

= E 2 



M l(r' r )M2(r' t )(-Q)-e 1 e2 Gm 1 m 2 -e 1 e 2 



(116) 



This approximation is always justified if the radius of the second intersection r" is 
enough large. Now, it is remarkable that whatever the value of r'* is, the first term 
in the square brackets in Eqn. (|116[) satisfies the inequality 



M 1 «)M 2 «)(-Q)-eie 2 Gm x m 2 - ei e 2 



(117) 



Comparing the expressions (|116jl . (|1 1T[1 and (j 1 12(1 it is possible to see that in the 
relativistic regime the shell ejection possibility is even greater than in the Newtonian 
case. Furthermore, it is worth noting that the presence of the charge do not change 
qualitatively the pure gravitational analysis, but just magnifies the ejection effect. 



7.1. Gm^m-2, — e^e^ < case 

Let us consider also briefly the case in which the shells are equal-signed charged 
and the repulsion overcome the gravity attraction, i.e. Gm\ra2 — e\e 2 < 0. In this 
case the ejection can happen only after an odd number of intersections. 

E.g. after three intersections, from the previous formulas we have, in the New- 
tonian approximation: 



£'" =£% - (Gmim 2 - e Y e 2 ) [ ~r - -jj + -ttt 



Obviously this formula has a meaning only if 

1 1 1 



(118) 



(119) 



e Remember that — Q = 1 + o(l/c 2 ). 
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otherwise the ejection happens at the first intersection (and then there would not be 
other crossings, and no r", r"'), or never more; if Eqn. (|119p is true, then it means 
that the barycenter of the two shells is falling into the center singularity. 



8. Conclusions 

We have found the energy exchange between two charged crossing shells (formula 
(|90| ). Then we have studied special cases of physical interest in which the formulas 
simplify: the non relativistic case, the massless shells, the test shell, and finally 
the ejection mechanism in a semi-Newtonian regime: we found that the ejection 
mechanism is more efficient in the charged case than in the neutral one if the 
charges have opposite sign (because the energy transfer is larger due to the Coulomb 
interaction). 
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